Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family.

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

Genetic improvement of Egyptian cotton (Gossypium barbadense L.) for high yield and fiber quality properties under semi arid conditions.

Between 2016 and 2018, the Agriculture Research Center's Sakha Agriculture Research Station conducted two rounds of pedigree selection on a segregating population of cotton (Gossypium barbadense L.) using the F2, F3, and F4 generations resulting from crossing Giza 94 and Suvin. In 2016, the top 5% of plants from the F2 population were selected based on specific criteria. The superior families from the F3 generation were then selected to produce the F4 families in 2017, which were grown in the 2018 summer season in single plant progeny rows and bulk experiments with a randomized complete block design of three replications. Over time, most traits showed increased mean values in the population, with the F2 generation having higher Genotypic Coefficient of Variance (GCV) and Phenotypic Coefficient of Variance (PCV) values compared to the succeeding generations for the studied traits. The magnitude of GCV and PCV in the F3 and F4 generations was similar, indicating that genotype had played a greater role than the environment. Moreover, the mean values of heritability in the broad sense increased from generation to generation. Selection criteria I2, I4, and I5 were effective in improving most of the yield and its component traits, while selection criterion I1 was efficient in improving earliness traits. Most of the yield and its component traits showed a positive and significant correlation with each other, highlighting their importance in cotton yield. This suggests that selecting to improveone or more of these traits would improve the others. Families number 9, 13, 19, 20, and 21 were the best genotypes for relevant yield characters, surpassing the better parent, check variety, and giving the best values for most characters. Therefore, the breeder could continue to use these families in further generations as breeding genotypes to develop varieties with high yields and its components.

Analysis of endophytic bacterial diversity in seeds of different genotypes of cotton and the suppression of Verticillium wilt pathogen infection by a synthetic microbial community.

BACKGROUND: In agricultural production, fungal diseases significantly impact the yield and quality of cotton (Gossypium spp.) with Verticillium wilt posing a particularly severe threat. RESULTS: This study is focused on investigating the effectiveness of endophytic microbial communities present in the seeds of disease-resistant cotton genotypes in the control of cotton Verticillium wilt. The technique of 16S ribosomal RNA (16S rRNA) amplicon sequencing identified a significant enrichment of the Bacillus genus in the resistant genotype Xinluzao 78, which differed from the endophytic bacterial community structure in the susceptible genotype Xinluzao 63. Specific enriched strains were isolated and screened from the seeds of Xinluzao 78 to further explore the biological functions of seed endophytes. A synthetic microbial community (SynCom) was constructed using the broken-rod model, and seeds of the susceptible genotype Xinluzao 63 in this community that had been soaked with the SynCom were found to significantly control the occurrence of Verticillium wilt and regulate the growth of cotton plants. Antibiotic screening techniques were used to preliminarily identify the colonization of strains in the community. These techniques revealed that the strains can colonize plant tissues and occupy ecological niches in cotton tissues through a priority effect, which prevents infection by pathogens. CONCLUSION: This study highlights the key role of seed endophytes in driving plant disease defense and provides a theoretical basis for the future application of SynComs in agriculture.

Characterization and expression analysis of GLABRA3 (GL3) genes in cotton: insights into trichome development and hormonal regulation.

BACKGROUND: GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) genes encode a typical helix-loop-helix (bHLH) transcription factors that primarily regulate trichome branching and root hair development, DNA endoreduplication, trichoblast size, and stomatal formation. The functions of GL3 genes in cotton crop have been poorly characterized. In this study, we performed comprehensive genome-wide scans for GL3 and EGL3 homologs to enhance our comprehension of their potential roles in trichome and fiber development in cotton crop. METHODS AND RESULTS: Our findings paraded that Gossypium hirsutum and G. barbadense have 6 GL3s each, unevenly distributed on 4 chromosomes whereas, G. arboreum, and G. raimondii have 3 GL3s each, unevenly distributed on 2 chromosomes. Gh_A08G2088 and Gb_A09G2187, despite having the same bHLH domain as the other GL3 genes, were excluded due to remarkable short sequences and limited number of motifs, indicating a lack of potential functional activity. The phylogenetic analysis categorized remaining 16 GL3s into three subfamilies (Group I-III) closely related to A. thaliana. The 16 GL3s have complete bHLH domain, encompassing 590-631 amino acids, with molecular weights (MWs) ranging from 65.92 to 71.36 kDa. Within each subfamily GL3s depicted shared similar gene structures and motifs, indicating conserved characteristics within respective groups. Promoter region analysis revealed 27 cis-acting elements, these elements were responsive to salicylic acid, abscisic acid (ABA), methyl jasmonate (MeJA), and gibberellin. The expression of GL3 genes was analyzed across 12 tissues in both G. barbadense and G. hirsutum using the publicly available RNA-seq data. Among GL3s, Gb_D11G0219, Gb_D11G0214, and Gb_D08G2182, were identified as relatively highly expressed across different tissues, consequently selected for hormone treatment and expression validation in G. barbadense. RT-qPCR results demonstrated significant alterations in the expression levels of Gb_D11G0219 and Gb_D11G0214 following MeJA, GA, and ABA treatment. Subcellular localization prediction revealed that most GL3 proteins were predominantly expressed in the nucleus, while a few were localized in the cytoplasm and chloroplasts. CONCLUSIONS: In summary, this study lays the foundation for subsequent functional validation of GL3 genes by identifying hormonal regulation patterns and probable sites of action in cotton trichome formation and fiber development. The results stipulate a rationale to elucidate the roles and regulatory mechanisms of GL3 genes in the intricate process of cotton fibre and trichome development.

Fingerprint Finder: Identifying Genomic Fingerprint Sites in Cotton Cohorts for Genetic Analysis and Breeding Advancement.

Genomic data in Gossypium provide numerous data resources for the cotton genomics community. However, to fill the gap between genomic analysis and breeding field work, detecting the featured genomic items of a subset cohort is essential for geneticists. We developed FPFinder v1.0 software to identify a subset of the cohort's fingerprint genomic sites. The FPFinder was developed based on the term frequency-inverse document frequency algorithm. With the short-read sequencing of an elite cotton pedigree, we identified 453 pedigree fingerprint genomic sites and found that these pedigree-featured sites had a role in cotton development. In addition, we applied FPFinder to evaluate the geographical bias of fiber-length-related genomic sites from a modern cotton cohort consisting of 410 accessions. Enriching elite sites in cultivars from the Yangtze River region resulted in the longer fiber length of Yangze River-sourced accessions. Apart from characterizing functional sites, we also identified 12,536 region-specific genomic sites. Combining the transcriptome data of multiple tissues and samples under various abiotic stresses, we found that several region-specific sites contributed to environmental adaptation. In this research, FPFinder revealed the role of the cotton pedigree fingerprint and region-specific sites in cotton development and environmental adaptation, respectively. The FPFinder can be applied broadly in other crops and contribute to genetic breeding in the future.

Integrated analysis of the transcriptome and metabolome reveals the molecular mechanism regulating cotton boll abscission under low light intensity.

BACKGROUND: Cotton boll shedding is one of the main factors adversely affecting the cotton yield. During the cotton plant growth period, low light conditions can cause cotton bolls to fall off prematurely. In this study, we clarified the regulatory effects of low light intensity on cotton boll abscission by comprehensively analyzing the transcriptome and metabolome. RESULTS: When the fruiting branch leaves were shaded after pollination, all of the cotton bolls fell off within 5 days. Additionally, H2O2 accumulated during the formation of the abscission zone. Moreover, 10,172 differentially expressed genes (DEGs) and 81 differentially accumulated metabolites (DAMs) were identified. A KEGG pathway enrichment analysis revealed that the identified DEGs and DAMs were associated with plant hormone signal transduction and flavonoid biosynthesis pathways. The results of the transcriptome analysis suggested that the expression of ethylene (ETH) and abscisic acid (ABA) signaling-related genes was induced, which was in contrast to the decrease in the expression of most of the IAA signaling-related genes. A combined transcriptomics and metabolomics analysis revealed that flavonoids may help regulate plant organ abscission. A weighted gene co-expression network analysis detected two gene modules significantly related to abscission. The genes in these modules were mainly related to exosome, flavonoid biosynthesis, ubiquitin-mediated proteolysis, plant hormone signal transduction, photosynthesis, and cytoskeleton proteins. Furthermore, TIP1;1, UGT71C4, KMD3, TRFL6, REV, and FRA1 were identified as the hub genes in these two modules. CONCLUSIONS: In this study, we elucidated the mechanisms underlying cotton boll abscission induced by shading on the basis of comprehensive transcriptomics and metabolomics analyses of the boll abscission process. The study findings have clarified the molecular basis of cotton boll abscission under low light intensity, and suggested that H2O2, phytohormone, and flavonoid have the potential to affect the shedding process of cotton bolls under low light stress.

A genome-wide association study for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in diploid cotton (Gossypium arboreum) and resistance transfer to tetraploid Gossypium hirsutum.

Fusarium wilt, caused by the soilborne fungus Fusarium oxysporum f. sp. vasinfectum (FOV), is a devastating disease affecting cotton (Gossypium spp.) worldwide. Understanding the genetic basis of resistance in diploid cotton and successfully transferring the resistance to tetraploid Upland cotton (G. hirsutum) are crucial for developing resistant cotton cultivars. Although numerous studies have been conducted to investigate the genetic basis of Fusarium wilt in tetraploid cotton, little research has been conducted on diploid species. In this study, an association mapping panel consisting of 246 accessions of G. arboreum, was used to identify chromosomal regions for FOV race 4 (FOV4) resistance based on foliar disease severity ratings in four greenhouse tests. Through a genome-wide association study (GWAS) based on 7,009 single nucleotide polymorphic (SNP) markers, 24 FOV4 resistance QTLs, including three major QTLs on chromosomes A04, A06, and A11, were detected. A validation panel consisting of 97 diploid cotton accessions was employed, confirming the presence of several QTLs. Evaluation of an introgressed BC2F7 population derived from G. hirsutum/G. aridum/G. arboreum showed significant differences in disease incidence and mortality rate, as compared to susceptible and resistant controls, suggesting that the resistance in G. arboreum and/or G. aridum was transferred into Upland cotton for the first time. The identification of novel major resistance QTLs, along with the transfer of resistance from the diploid species, expands our understanding of the genomic regions involved in conferring resistance to FOV4 and contributes to the development of resilient Upland cotton cultivars.

A systematic analysis of ARM genes revealed that GhARM144 regulates the resistance against Verticillium dahliae via interaction with GhOSM34.

Proteins of the armadillo repeat gene family play important roles in plant pathogen response. Here, 169 armadillo (ARM) genes were identified in upland cotton (Gossypium hirsutum). Phylogenetic analysis grouped these into 11 subfamilies, with conserved protein structures within each subfamily. The results signify that the expansion of the gene family occurred via whole genome duplication and dispersed duplication. Expression profiling and network analysis suggest that GhARM144 may regulate cotton resistance to Verticillium dahliae. GhARM144 was upregulated in roots by V. dahliae infection or salicylic acid treatment. This upregulation indicates a negative regulatory role of GhARM144' in the cotton immune responses, potentially by manipulating salicylic acid biosynthesis. Protein interaction studies found that GhARM144 associates with an osmotin-like protein, GhOSM34, at the plasma membrane. Silencing GhOSM34 reduced the resistance to V. dahliae, suggesting it may play a positive regulatory role. The results demonstrate that GhARM144 modulates cotton immunity through interaction with GhOSM34 and salicylic acid signalling. Further study of these proteins may yield insights into disease resistance mechanisms in cotton and other plants.

Genome-wide identification of the key kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L.).

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (G. hirsutum L.). Results showed that 159 kinesin genes, including 15 genes of the kinesin-13 gene subfamily, were identified in upland cotton; of which 157 kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 kinesin genes in upland cotton, including 10 kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 kinesin genes were significantly associated with three fiber traits, among which a kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

AmCBF1 activates the expression of GhClpR1 to mediate dark-green leaves in cotton (Gossypium hirsutum).

KEY MESSAGE: The transcription factor AmCBF1 deepens the leaf colour of transgenic cotton by binding to the promoter of the chloroplast development-related protein GhClpR1 to promote photosynthesis. The ATP-dependent caseinolytic protease (Clp protease) family plays a crucial role within chloroplasts, comprising several Clp proteins to maintain chloroplast homeostasis. At present, research on Clp proteins mainly focuses on Arabidopsis, leaving its function in other plants, particularly in crops, less explored. In this study, we overexpressed AmCBF1 from Ammopiptanthus mongolicus (A. mongolicus) in wild type (R15), and found a significant darkening of leaf colour in transgenic plants (L28 and L30). RNA-seq analysis showed an enrichment of pathways associated with photosynthesis. Subsequent screening of differentially expressed genes revealed a significant up-regulation of GhClpR1, a gene linked to chloroplast development, in the transgenic strain. In addition, GhClpR1 was consistently expressed in upland cotton, with the highest expression observed in leaves. Subcellular localization analysis revealed that the protein encoded by GhClpR1 was located in chloroplasts. Yeast one hybrid and dual luciferase experiments showed that the AmCBF1 transcription factor positively regulates the expression of GhClpR1. VIGs-mediated silencing of GhClpR1 led to a significant yellowing phenotype in the leaves. This was accompanied by a reduction in chlorophyll content, and microscopic examination of chloroplast ultrastructure revealed severe developmental impairment. Finally, yeast two-hybrid assays showed that GhClpR1 interacts with the Clp protease complex accessory protein GhClpT2. Our study provides a foundation for studying the function of the Clp protease complex and a new strategy for cultivating high-light-efficiency cotton resources.

An insight into the gene expression evolution in Gossypium species based on the leaf transcriptomes.

BACKGROUND: Gene expression pattern is associated with biological phenotype and is widely used in exploring gene functions. Its evolution is also crucial in understanding species speciation and divergence. The genus Gossypium is a bona fide model for studying plant evolution and polyploidization. However, the evolution of gene expression during cotton species divergence has yet to be extensively discussed. RESULTS: Based on the seedling leaf transcriptomes, this work analyzed the transcriptomic content and expression patterns across eight cotton species, including six diploids and two natural tetraploids. Our findings indicate that, while the biological function of these cotton transcriptomes remains largely conserved, there has been significant variation in transcriptomic content during species divergence. Furthermore, we conducted a comprehensive analysis of expression distances across cotton species. This analysis lends further support to the use of G. arboreum as a substitute for the A-genome donor of natural cotton polyploids. Moreover, our research highlights the evolution of stress-responsive pathways, including hormone signaling, fatty acid degradation, and flavonoid biosynthesis. These processes appear to have evolved under lower selection pressures, presumably reflecting their critical role in the adaptations of the studied cotton species to diverse environments. CONCLUSIONS: In summary, this study provided insights into the gene expression variation within the genus Gossypium and identified essential genes/pathways whose expression evolution was closely associated with the evolution of cotton species. Furthermore, the method of characterizing genes and pathways under unexpected high or slow selection pressure can also serve as a new strategy for gene function exploration.

GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, play a regulatory role in cotton fiber secondary cell wall biosynthesis.

KEY MESSAGE: GhHB14_D10 and GhREV_D5 regulated secondary cell wall formation and played an important role in fiber development. Cotton serves as an important source of natural fiber, and the biosynthesis of the secondary cell wall plays a pivotal role in determining cotton fiber quality. Nevertheless, the intricacies of this mechanism in cotton fiber remain insufficiently elucidated. This study investigates the functional roles of GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, in secondary cell wall biosynthesis in cotton fibers. Both GhHB14_D10 and GhREV_D5 were found to be localized in the nucleus with transcriptional activation activity. Ectopic overexpression of GhHB14_D10 and GhREV_D5 in Arabidopsis resulted in changed xylem differentiation, secondary cell wall deposition, and expression of genes related to the secondary cell wall. Silencing of GhHB14_D10 and GhREV_D5 in cotton led to enhanced fiber length, reduced cell wall thickness, cellulose contents and expression of secondary cell wall-related genes. Moreover, GhHB14_D10's direct interaction with GhREV_D5, and transcriptional regulation of cellulose biosynthesis genes GhCesA4-4 and GhCesA7-2 revealed their collaborative roles in secondary cell wall during cotton fiber development. Overall, these results shed light on the roles of GhHB14_D10 and GhREV_D5 in secondary cell wall biosynthesis, offering a strategy for the genetic improvement of cotton fiber quality.

Melatonin-mediated regulation of autophagy is independent of ABA under drought stress in sensitive variety of Gossypium hirsutum L.

Autophagy is a highly conserved process that plays a crucial role in adaptation of plants to stress conditions. Melatonin and abscisic acid (ABA) share an antagonistic relationship; however, both are reported to elevate autophagy individually. Here, we report that melatonin alleviates drought stress effects like wilting and stunted growth in 18-day-old plants of drought-sensitive variety of cotton (Gossypium hirsutum L.) and improves the plant growth, chlorophyll content, photosynthetic efficiency, and sugar metabolism and transport. Melatonin priming increased the endogenous melatonin content (5.02-times) but decreased the ABA (2.63-times) by reducing NCED3 expression as compared to unprimed plants under drought. Also, elevated expression of ATG8c and ATG8f correlated with higher lipidated-ATG8 levels and modulation of RAPTOR1 suggesting a higher occurrence of autophagy and regulation of plant growth in primed stressed plants. Additionally, decreased TPS63 and increased TPP22 expression could have lowered the accumulation of trehalose-6-P (T6P) in primed stressed plants thus contributing to autophagy progression. Priming also enhanced the expression of MAPK6 and RAF18, and increased the transcript/protein levels of SnRK2.6 and KIN10, which is pointing towards melatonin's beneficial effect on autophagy under drought. Despite higher ABA content, elevated TPS63 and downregulated TPP22 could have hindered autophagy induction in unprimed stressed plants. Although fluridone treatment reduced the ABA content, the expression of SnRK2.6 and KIN10 remained unaltered in fluridone-treated and untreated primed plants indicating the ABA-independent expression. These results suggest that the melatonin-mediated activation of MAPK contributes to the ABA-independent activation of SnRK2, consequently, SnRK1 and autophagy under drought.

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.).

KEY MESSAGE: GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

The comparison of species diversity and abundance of insect natural enemies in the domesticated species of cotton using the yellow pan trap method.

India is the world's largest cotton producer and the only country that grows all four cultivated cotton species. There have been very few studies on the diversity and abundance of natural enemies of cotton insect pests in these cultivated cotton species. Therefore, the current study (2016-2018) was conducted to assess the diversity and abundance of natural enemies that cultivated cotton species harbour. Phule Dhanwantari, Suraj, Suvin, RCH-2, and DCH-32 were the five genotypes used in the study, each with a distinct genetic background. Using the adiv 2.0.1 and vegan R packages, we identified significant differences in natural enemies in terms of species diversity, richness, evenness, and abundance. Analysis of Similarity (ANOSIM) and Non-metric Multidimensional Scaling (NMDS) indicated substantial differences in the natural enemy community structure among the examined genotypes. A total of 17,279 natural enemies were collected and identified across genotypes from seven predatory families and five parasitoid families. The percentage share of these natural enemy families across genotypes and years, in descending order, is Coccinellidae (28.23%) < Tachinidae (19.23%) < Braconidae (12.68%) < Chrysopidae (11.65%) < Chalcididae (9.41%) < Aphelinidae (6.33%) < Pentatomidae (3.29%) < Ichneumonidae (2.37%) < Syrphidae (2.33%) < Vespidae (1.81%) < Asilidae (1.79%) < Geocoridae (0.89%). Coccinellidae, Tachinidae, Braconidae, Chrysopidae, Chalcididae, and Aphelinidae are the six major families that account for more than 85% of all recorded natural enemies. These six families have a higher percentage share in Phule Dhanwantary (90%) compared to the other genotypes. The conservation and better utilization of these natural enemies are crucial for the ecological and safe management of insect pests in the cotton ecosystem.

Rapidly mining candidate cotton drought resistance genes based on key indicators of drought resistance.

BACKGROUND: Focusing on key indicators of drought resistance is highly important for quickly mining candidate genes related to drought resistance in cotton. RESULTS: In the present study, drought resistance was identified in drought resistance-related RIL populations during the flowering and boll stages, and multiple traits were evaluated; these traits included three key indicators: plant height (PH), single boll weight (SBW) and transpiration rate (Tr). Based on these three key indicators, three groups of extreme mixing pools were constructed for BSA-seq. Based on the mapping interval of each trait, a total of 6.27 Mb QTL intervals were selected on chromosomes A13 (3.2 Mb), A10 (2.45 Mb) and A07 (0.62 Mb) as the focus of this study. Based on the annotation information and qRT‒PCR analysis, three key genes that may be involved in the drought stress response of cotton were screened: GhF6'H1, Gh3AT1 and GhPER55. qRT‒PCR analysis of parental and extreme germplasm materials revealed that the expression of these genes changed significantly under drought stress. Cotton VIGS experiments verified the important impact of key genes on cotton drought resistance. CONCLUSIONS: This study focused on the key indicators of drought resistance, laying the foundation for the rapid mining of drought-resistant candidate genes in cotton and providing genetic resources for directed molecular breeding of drought resistance in cotton.

An Oxidoreductase-like Protein is Required for Verticillium dahliae Infection and Participates in the Metabolism of Host Plant Defensive Compounds.

Verticillium dahliae, a notorious phytopathogenic fungus, is responsible for vascular wilt diseases in numerous crops. Uncovering the molecular mechanisms underlying pathogenicity is crucial for controlling V. dahliae. Herein, we characterized a putative oxidoreductase-like protein (VdOrlp) from V. dahliae that contains a functional signal peptide. While the expression of VdOrlp was low in artificial media, it significantly increased during host infection. Deletion of VdOrlp had minimal effects on the growth and development of V. dahliae but severely impaired its pathogenicity. Metabolomic analysis revealed significant changes in organic heterocyclic compounds and phenylpropane compounds in cotton plants infected with ΔVdOrlp and V991. Furthermore, VdOrlp expression was induced by lignin, and its deletion affected the metabolism of host lignin and phenolic acids. In conclusion, our results demonstrated that VdOrlp plays an important role in the metabolism of plant phenylpropyl lignin and organic heterocyclic compounds and is required for fungal pathogenicity in V. dahliae.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum.

KEY MESSAGE: GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Genome-Wide Identification of the GhANN Gene Family and Functional Validation of GhANN11 and GhANN4 under Abiotic Stress.

Annexins (ANNs) are a structurally conserved protein family present in almost all plants. In the present study, 27 GhANNs were identified in cotton and were unevenly distributed across 14 chromosomes. Transcriptome data and RT-qPCR results revealed that multiple GhANNs respond to at least two abiotic stresses. Similarly, the expression levels of GhANN4 and GhANN11 were significantly upregulated under heat, cold, and drought stress. Using virus-induced gene silencing (VIGS), functional characterization of GhANN4 and GhANN11 revealed that, compared with those of the controls, the leaf wilting of GhANN4-silenced plants was more obvious, and the activities of catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD) were lower under NaCl and PEG stress. Moreover, the expression of stress marker genes (GhCBL3, GhDREB2A, GhDREB2C, GhPP2C, GhRD20-2, GhCIPK6, GhNHX1, GhRD20-1, GhSOS1, GhSOS2 and GhSnRK2.6) was significantly downregulated in GhANN4-silenced plants after stress. Under cold stress, the growth of the GHANN11-silenced plants was significantly weaker than that of the control plants, and the activities of POD, SOD, and CAT were also lower. However, compared with those of the control, the elasticity and orthostatic activity of the GhANN11-silenced plants were greater; the POD, SOD, and CAT activities were higher; and the GhDREB2C, GhHSP, and GhSOS2 expression levels were greater under heat stress. These results suggest that different GhANN family members respond differently to different types of abiotic stress.

Isolation and Functional Characterization of a Constitutive Promoter in Upland Cotton (Gossypium hirsutum L.).

Multiple cis-acting elements are present in promoter sequences that play critical regulatory roles in gene transcription and expression. In this study, we isolated the cotton FDH (Fiddlehead) gene promoter (pGhFDH) using a real-time reverse transcription-PCR (qRT-PCR) expression analysis and performed a cis-acting elements prediction analysis. The plant expression vector pGhFDH::GUS was constructed using the Gateway approach and was used for the genetic transformation of Arabidopsis and upland cotton plants to obtain transgenic lines. Histochemical staining and a ß-glucuronidase (GUS) activity assay showed that the GUS protein was detected in the roots, stems, leaves, inflorescences, and pods of transgenic Arabidopsis thaliana lines. Notably, high GUS activity was observed in different tissues. In the transgenic lines, high GUS activity was detected in different tissues such as leaves, stalks, buds, petals, androecium, endosperm, and fibers, where the pGhFDH-driven GUS expression levels were 3-10-fold higher compared to those under the CaMV 35S promoter at 10-30 days post-anthesis (DPA) during fiber development. The results indicate that pGhFDH can be used as an endogenous constitutive promoter to drive the expression of target genes in various cotton tissues to facilitate functional genomic studies and accelerate cotton molecular breeding.